Characterizing the Hofstadter butterfly's outline with Chern numbers 
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In this work, we report original properties inherent to independent particles subjected to a magnetic field 
by emphasizing the existence of regular structures in the energy spectrum's outline. We show that this fractal 
curve, the well-known Hofstadter butterfly's outline, is associated to a specific sequence of Chern numbers 
that correspond to the quantized transverse conductivity. Indeed the topological invariant that characterizes the 
fundamental energy band depicts successive stairways as the magnetic flux varies. Moreover each stairway is 
shown to be labeled by another Chern number which measures the charge transported under displacement of the 
periodic potential. We put forward the universal character of these properties by comparing the results obtained 
for the square and the honeycomb geometries. 
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I. INTRODUCTION 

The energy spectrum associated to charged particles 
moving in a two-dimensional lattice and subjected to a 
high magnetic field has inspired numerous works, since the 
seminal papers of Hofstadter HI and Wannier [2]. When the 
magnetic flux penetrating the lattice is a rational number, 
namely $ = p/q where p, q are integers, the spectrum splits 
into q subbands. The representation of the spectrum as a 
function of the flux shows a recursive structure with clear 
self-similarities [3|, and adopts the shape of an intriguing 
insect: the so-called Hofstadter butterfly [1]. The infinitely 
many gaps which compose this surprising figure are known 
to follow a simple rule: each gap is labeled by two integers 
(t r ,s r ) which satisfy a Diophantine equation (21 El- In a 
fundamental work, Thouless et al. have emphasized that 
the Hofstadter butterfly plays a key role in the quantum 
Hall effect theory: when the Fermi energy of the system 
lies in a gap, the transverse conductivity of the system is 
quantized and is given by a xy = (e 2 /h) t r , where e is the 
particle's charge, h is Planck's constant and t r satisfies the 
aforementioned Diophantine equation ||6) . 
In this context, the Green-Kubo expression for the transverse 
conductivity has an elegant topological interpretation: when 
the Fermi energy of the system lies in a gap, the transverse 
conductivity is given by a sum of Chern numbers. The latter 
are associated with the filled bands that are situated below 
the gap (4). From a mathematical point of view, the Chern 
number is an integral invariant characterizing the topology of 
a fibre bundle which is defined on each energy band j4l [7). 
Moreover, these topological invariants are known to be the 
only quantized quantities which can be associated with the 
energy bands |8|. The Diophantine equation, which describes 
any 2D electronic systems subjected to a magnetic field J5), 
is highly connected to topology since both integers (t r ,s r ) 
can be interpreted as Chern numbers (9][T0]. As pointed out 
by MacDonald and Kunz, the Chern number s r measures the 
charge transported when the periodic potential is adiabatically 
displaced J3]|9l- 

Recent advances in cold atoms physics ifTTl [T2l allow for 



the experimental exploration of the Hofstadter butterfly using 
atoms trapped in optical lattices [ 13 1. Different arrangements, 
which indeed mimic the presence of an artificial "magnetic" 
field in the dynamics of neutral particles lfT3l[T4l[T5l[T6l[T7l . 
are nowadays studied in laboratories ifTSl . With such se- 
tups one expect to explore new features in the fields of 
vortex physics lfT2l [T9l l20l |2D and quantum Hall systems 
El |22l EH 131 125] |26). In particular, we have suggested 
that an integer quantum Hall-like effect for neutral fermionic 
particles should be observed in optical lattices 11221 . 
Recently, we have also shown that fermionic atoms trapped 
in 2D optical lattices and subjected to an artificial "magnetic" 
field should undergo a Mott metal-insulator transition and that 
the phase boundary depicts the Hofstadter butterfly's outline 
ETll . Although the Mott-insulator phase transition occurs 
in the system when the interaction between the particles is 
taken into account, the phase boundary only depends on the 
underlying single-particle physics. A similar result has been 
obtained by Oktel et al. Ctl|29] and by Goldbaum et al. J30| 
in the context of the bosonic superfluid-insulator transition. 
Motivated by the important role played by the Hofstadter 
butterfly's outline in this theoretical framework, and by 
the recent experimental advances in the field of ultracold 
atoms, we investigate intrinsic properties associated to this 
fractal and intriguing curve. We point out that the Hofstadter 
butterfly's outline is also known to represent the phase 
boundary for the normal-superconducting phase transition in 
superconducting networks OTl l32l [33l . 



In this paper, we characterize the Hofstadter butterfly's out- 
line by studying the Chern numbers associated to the funda- 
mental energy band of the system as a function of the mag- 
netic flux. Physically, we evaluate the system's transverse 
conductivity when the Fermi energy lies in the first gap of 
the energy spectrum. Under this condition, we show that the 
measure of the transverse conductivity exhibits a particular 
sequence of plateaus as a function of the magnetic flux. We 
study such structures and emphasize their universal character 
by comparing the results obtained for the square and for the 
honeycomb geometries. 
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II. THE SQUARE LATTICE 

We first consider the case of a fermionic gas trapped in a 
2D optical square lattice and subjected to an artificial "mag- 
netic" field lfl2l IT~3l . We assume that the optical potential 
created by the lasers is sufficiently strong in order to apply 
a tight-binding approximation and we investigate this system 
in the non-interacting limit, which can be reached at low den- 
sities. The magnetic field B = B l z , which is characterized 
by the number of effective magnetic flux quanta per unit cell 

is supposed to be induced by rotation lfT9ll25l . or by com- 
bining laser-assisted tunneling and lattice acceleration meth- 
ods iTPJl IT7l [181 or by immersion of the lattice into a rotating 
Bose-Einstein condensate [16]. For a rotating system, this pa- 
rameter is given by $ = 2Ma 2 il/h, where f2 is the angular 
velocity, M is the particle's mass, a is the lattice constant and 
h = 2ttH is Planck's constant ll25ll . We can treat this problem 
by considering the Landau gauge, A = (0, Bx, 0), for which 
the many-body Hamiltonian reads H = t J2 m n ^m,n> with 

^m,n — ^m,n^^ — ^ ^m.n^^^n—l ~t~ b-.C, (1) 

where a m , n (a^ n ) is the fermionic annihilation (creation) op- 
erator on the lattice site (m, n) and t is the nearest-neighbor 
tunneling amplitude. In this work, the effective magnetic flux 
$ is supposed to vary between [0, 1], which is achieved for 
w 1 kHz. In the following, we set the lattice constant a and 
the particle's mass M to unity and work in units where h = 1 
except otherwise stated. For rational fluxes, $ = p/q where 
p and q are integers, the single-particle Schrodinger equation 
yields the well-known Harper equation 0] 

e lK -0m+i + e~ lfcx Vm-i + 2cos(27r$m-fc y )^ m = — ip m , 

(2) 

where ip m is a g-periodic wave function, fc is the wave vec- 
tor and E is the single particle energy. The wave vector be- 
longs to the magnetic Brillouin zone, a two-torus defined as 
k x £ [0, ~] and k y £ [0, 2ir]. The energy spectrum asso- 
ciated to Eq. |2]) has a band structure, composed of q sub- 
bands, which has been extensively studied in the literature 
lfTll2l[3l l34ll35ll36ll : the representation of the energy as a func- 
tion of the flux leads to the Hofstadter butterfly (Fig jTJ. This 
striking fractal figure, which illustrates the infinitely many 
gaps of the spectrum, exhibits a recursive structure and is de- 
scribed through a simple rule Q: for $ = p/q, the r th gap of 
the spectrum is labeled by two integers (i r , s r ), which satisfy 
a Diophantine equation Q [5) 

r = pt r + qs r . (3) 

In the square lattice case, the condition \t r \ < q/2 determines 
the solution unambiguously • 



In this work, we focus our attention on the outline of the 
butterfly (red curve in Fig{T|, which is known to play a key 
role in the field of quantum phase transitions 1 27 28 29 30 




E 



Figure 1: (Color online) Hofstadter butterfly $ = &(E): single- 
particle spectrum for a 2D square lattice subjected to a magnetic 
field. The butterfly's outline, which corresponds to the fundamental 
energy band Ei(k x , k y ), is colored in red. The energy is expressed 
in units of the hopping parameter t. 

|3T1 l32l l33l . Our aim is to investigate the structures charac- 
terizing this non-trivial curve, which is formed by the fun- 
damental band Ei(k x , k y ), where H\%jji) = Ei\ip±), in the 
range ^ — p/q £ [0, 1]. As q increases, the butterfly's out- 
line doesn't seem to smoothen because of its fractal nature [ 1] . 
In order to maintain the system along this irregular curve, we 
suppose that the Fermi energy lies inside the first gap of the 
single particle spectrum for all $ £ [0,1]. Under this assump- 
tion, we are able to compute the analogue of the transverse 
conductivity for neutral currents l22l . which is associated to 
the energy band Ei(k x , k y ) for a given value of the flux 
In this context, the transverse conductivity can be computed 
with Kubo's formula and is expressed as l6l l36ll 

cr = 1 1 [ (< 9 ^| <9 ^ _ \<ttl\ (4) 

xy h 2wi Jj2 dk x dk y dk y dk x ' 

where the fundamental state \ip±) alone contributes. 

The quantization of this quantity follows from the topolog- 
ical interpretation of Eq. Q: the transverse conductivity is 
related to the topologically invariant Chem number C$ 071 . 
an integer defined as 

S 27T 

= — [ (^i - l^tl 1^-) (5) 

2-7T Jf2 dk x dk y dky dk x ' 

where T is the so-called Berry's curvature associated to the 
band Ei(k x , k y ). According to Refs. |6, 36], one finds that 
C s = —ti, and therefore this invariant integer satisfies the 
Diophantine equation Eq.([3| with r = 1, 

C s = $- 1 s 1 --, (6) 
P 
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Figure 2: Chem numbers as a function of the magnetic flux, \ Cs\ = 
|C s |($),for$ = 2 withq < 77. 



with the condition \Cs\ < <?/2. 

One can solve Eq.|6]) for all <f> = p/q S [0, 1], in order 
to obtain the many Chern numbers associated to the Hofs- 
tadter butterfly's outline. Technically one fixes a high value 
for the denominator q and computes the Chern number Cs for 
p = 1, 2, q, such that p and g are mutually primes. 
The illustration of these integers as a function of the effec- 
tive magnetic flux is quite surprising. All along the irregular 
outline E\ = Ex{&), the Chern numbers computed for the 
various <f> follow a very regular law: the representation of the 
Chern numbers as a function of the flux Cs = Cs(&) depicts 
plateau sequences, adopting the shape of successive stairways. 
In Fig|2j we show this structure in a compact way, by plotting 
\Cs\ as a function of the flux We note that this figure is 
symmetric with respect to $ = 0.5. 

The numerical result illustrated in Fig. [2] is already inter- 
esting since it underlines the complexity of the Hofstadter 
butterfly's outline. A topological argument stipulates that 
if the fundamental band Ei(k x ,k y ) is well separated from 
the first excited state E 2 (k x , k y ), for all $ G [0, 1], then the 
Chern number Cs should remain constant |8|. However, this 
is not the case since we have shown that this topological 
number takes many different values as the flux varies. For q 
fixed and p = 1, ...q, we obtain that the number of different 
values Nq increases linearly with respect to the denominator 
q. We eventually note that the Chern number's magnitude 
\Cs\ takes the value of all the natural numbers up to Nc, 
namely \C S \ = 1, 2, N c , and that N c ~ q/2. Since q 
may take an arbitrarily high value, our calculations show that 
the first gap closes infinitely many times as the flux is varied. 
These numerous gap closings are illustarted in Fig. [3] where 
the two lowest energy bands are plotted as a function of the 
flux <f>. 



0.5r 



0.4- 



0.3 

© 



02 



0.1- 



°4 3 2 1 1 

E 

Figure 3: Two lowest energy bands, E\ (blue lines) and E2 (red 
lines), as a function of the flux $. The energy is expressed in units 
of the hopping parameter t. 



Furthermore, we find that a strong structure is hidden be- 
hind the intriguing Fig|2] We show that the other Chern num- 
ber which satisfies the Diophantine equation with t\ = —Cs, 

s 1 = 1 +$C s , (7) 

q 

classifies the many points \Cs\) plotted in Fig(2] in a 
rigorous manner. We illustrate this result in Fig|4j in which 
we connect all the points (<T>, \Cs\) that are associated to the 
same number |si| with a colored line. As \si\ increases, the 
color changes progressively from red to purple. It is very 
clear from this colored figure that all the plateaus belong 
to a specific stairway, labeled by |si|, since no crossing 
between the lines is observed. Moreover, we point out that 
the successive stairways correspond to increasing values of 
|si|. We note that the integer t\ is related to the quantized 
transverse conductivity, whereas the number s\ actually 
measures the charge transported when the periodic potential 
is adiabatically displaced (3j[9], 

In summary, the first gap of the spectrum closes infinitely 
many times as the magnetic flux is varied. A signature of these 
gap closings is given by the Chern number C$, which is asso- 
ciated to the first energy band, and which takes infinitely many 
different values as a function of the flux. Surprisingly, these 
many values give rise to a structure consisting of successive 
stairways which are characterized by the number s\. 



III. THE HONEYCOMB LATTICE 

In this section, we consider the case of a 2D honeycomb lat- 
tice subjected to a magnetic field. The honeycomb has a bipar- 
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Figure 4: (Color online) Chem numbers (black dots) as a function 
of the magnetic flux, \C S \ = jC s j($), for $ = | with q < 157. 
The dots are connected by colored lines according to the topological 
invariant | si | to which they are associated through the Diophantine 
equation Eq.(6]l. Note that the Chern numbers are single- valued. 



Figure 5: (Color online) Butterfly spectrum $ = &(E) for the 2D 
honeycomb lattice subjected to a magnetic field. The butterfly's out- 
line is colored in red. The energy is expressed in units of the hopping 
parameter t. 



tite structure and it is common to define two fermion operators 
a) A {r) and a) B (r), where r = mei + ne 2 (see for exemple 
Ref. 11381 ). The unit vectors are chosen as e.\ — (3/2, V3/2) 
and e-2 = (0,\/3). In the tight-binding approximation, the 
many-body Hamiltonian reads H = t J2 r T~L( r )> where 

H(r) = a\(r)a B (r) + e a ^ m a\{r)a B {r - e 2 ) 

+ a A (r + ei)a B (r) +h.c, (8) 

and $ = p/q is the effective magnetic flux quanta per unit 
cell. The single-particle Schrodinger equation associated to 
Eq.([8]) yields 

^1 + e +t2 ^ m - tk y^ A (m) + e lk ^ A {m + 1) = jip B (m), 

(l + e- i2 ^ m+ik ^ip B (m)+e- ik *i> B (m-l) = -^(m), 

(9) 

where ifjAi'm) and ip B (m) are q-periodic wave functions, k 
is the wave vector and E is the single-particle energy. The 
wave vector belongs to the magnetic Brillouin zone defined 
as k x £ [0, and k y € [0,27r]. The spectrum is depicted 
in Fig j5] as a function of the effective magnetic flux <I>, and 
illustrates a modified version of the Hofstadter butterfly. The 
band structure associated to Eq.|9) has been extensively stud- 
ied in order to investigate the very rich physics of graphene 
ll38l [39l l40l . The Chern numbers associated to the energy 
bands have been computed numerically by Hatsugai et al. 
and give rise to the anomalous quantum Hall effect: around 
E = 0, the transverse conductivity evolves by steps according 
to a xy — ±(2N + l)e 2 /h where N is an integer ||38l . 



The outline of the honeycomb butterfly (red curve in 
Pig |5j» is highly irregular and differs from the square lattice 
case in regards to its general shape. In order to characterize 
this other irregular curve, and compare it to the Hofstadter 
butterfly's outline, one has to evaluate the Chern numbers 
Ch associated to the fundamental band Ei(k x ,k y ) for the 
honeycomb lattice case. Unfortunately, contrary to the square 
lattice case, no Diophantine equation is known to describe 
the entire honeycomb butterfly [38 1: the gaps are labeled by 
two integers (t r , s r ) which satisfy Eq. (|3), but the condition 
\t r \ < q/2 is not always satisfied. Thus the Chern numbers 
cannot be unambiguously determined on the basis of the 
Diophantine equation. In the lack of such an equation, one 
has to compute the Chern numbers numerically. This can be 
easily achieved thanks to an efficient method developed by 
Fukui et al. ATI . 

The numerical results are shown in Figj6] A new sequence 
of plateaus is observed in the honeycomb case, for which 
the transverse conductivity evolves by steps according to 
a xy = ±(2N + l)e /h, where N is an integer. This is 
already surprising, because the fundamental band Ei(k x , k y ) 
corresponds to energies which are far from E = 0, where 
the anomalous behavior of double steps is expected l38l . 
From this observation, it seems that the main signature of 
the anomalous quantum Hall effect is already contained in 
the edge of the butterfly: if the Fermi energy of the system 
lies in the first gap while varying one should observe this 
fascinating property proper to graphene. From a topological 
point of view, this particular sequence of Chern numbers 
suggests that the honeycomb butterfly's outline radically 
differs from the Hofstadter butterfly's outline. 

However, it is worth noticing that the general structure de- 
picted in Figj6]is remarkably similar to the results shown in 
Fig|2]and Fig]4]in the context of the square lattice case. The 
universality arising from these figures indicates that the Dio- 



5 



14 

12 
10 
8 
6 
4 
2 




0.4 0.6 



Figure 6: (Color online) Chem numbers (black dots) as a function 
of the magnetic flux, \C H \ = \C H \(®), for $ = | with q < 37. 
The dots are connected by colored lines according to the topological 
invariant jsij to which they are associated through the Diophantine 
equation Eq.{3}. Note that the Chern numbers are single- valued. 



phantine equation also plays a key role in the honeycomb sys- 
tem. 

Although the condition \Cr\ < q/2 is not fulfilled, we have 
verified that all the solutions Ch satisfy the Diophantine 
equation Eq. Q, and that each stairway is again characterized 
by the other number | si | . 

In order to emphasize the universality suggested by Fig|2]and 
Figj6] we have connected with a colored line the solutions 
Cjf ($) according to their associated number |si| (see Figj6|. 
As in the square lattice case, the successive stairways are in- 
deed characterized by increasing values of the Chern number 

I'll- 



IV. CONCLUSION 



In this work, we have put forward the existence of strong 
and regular structures associated to the highly irregular 
Hofstadter butterfly's outline. These structures, which arise 
in square and honeycomb lattices, have a universal character 
and are related to the underlying topology of the system's 
fundamental energy band. The topology of the system is 
designated through the Chem number which gives its value 
to the quantized transverse conductivity. The correspon- 
dence between the results obtained for the square and the 
honeycomb lattices has been confirmed through numerical 
computations of the Chern numbers in the honeycomb 
case. We have verified that these solutions indeed satisfy 
the general Diophantine equation. It has been shown, for 
both geometries, that when the Fermi energy remains in 
the first gap, the transverse conductivity is highly irregular 
but evolves on stairways labeled by the other number that 
satisfies the Diophantine equation. We believe that these 
structures might play a role in the Mott-insulator transitions 
observed in rotating atomic systems, where single-particle 
properties are known to be dominant. We eventually point out 
that the properties emphasized in this work are not restricted 
to the field of cold atoms physics and could also be found in 
superconducting networks or 2D electronic systems. 
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